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А.Т.РЫБАК 


ОБЪЕМНАЯ ЖЕСТКОСТЬ И ЕЁ ВЛИЯНИЕ НА ДИНАМИКУ 
ГИДРОМЕХАНИЧЕСКОЙ СИСТЕМЫ 


Предлагается понятие объёмной жёсткости элементов гидравлической системы. 
Приводятся полученные аналитическим путём формулы для определения приве- 
дённых объёмных жёсткостей цилиндрической трубы, поршневой и штоковой поло- 
стей гидроцилиндра, плунжерного гидроцилиндра, а также формула для определе- 
ния приведённой объёмной жёсткости произвольного числа гидравлических эле- 
ментов, включённых в гидросистему параллельно. Обосновывается вывод о воз- 
можности гидромеханического резонанса, вызываемого переменной объёмной 
жёсткостью гидравлической системы. 

Ключевые слова: гидросистема, нестационарный режим, объёмная жёсткость, 
гидромеханический резонанс. 


Постановка задачи. Подавляющее большинство современных исследова- 
ний, посвящённых изучению динамики гидромеханических систем (ГМС), 
основаны на использовании теории автоматического регулирования и 
управления. Математическая модель гидравлической системы в этом слу- 
чае составляется с использованием уравнения неразрывности потока, при 
этом учёт сжимаемости жидкости и деформации трубопроводов осуществ- 
ляется путём введения в уравнение некоего расхода жидкости, величина 
которого зависит от приведённого модуля упругости трубопровода [1-4]. 
Однако в последние годы всё более широкое распространение получают 
методы расчёта силового гидравлического привода, учитывающие наруше- 
ние уравнения неразрывности во время его работы в неустановившемся 
режиме. В этом случае в математической модели взамен приведённого мо- 
дуля упругости системы удобнее использовать её приведённую объёмную 
жёсткость. 

Объёмная жёсткость жидкости. Изменение давления жидкости в раз- 
личных точках гидравлической системы можно проследить, обратившись к 
обобщённому закону Гука [5], из которого следует, что 


ж 


Е 
ар=-Е ы: ‚ или обозначив С =—^, (1) 
% У 
получим ар =-—С„АУ», (2) 


где И - начальный объём жидкости; @И, — приращение объёма жидкости, 
соответствующее приращению давления на 4р; Ё„ = Е. — объём- 
ж 
ный модуль упругости жидкости; В» — коэффициент объемного сжа- 
тия жидкости; Ск — объемная жёсткость жидкости, которая представ- 
ляет собой приращение давления жидкости, соответствующее еди- 
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ничному приращению первоначального объёма, или иначе, прираще- 
ние давления, необходимое для изменения первоначального объёма 
жидкости И на единицу. Знак «минус» указывает на то, что положи- 
тельному приращению давления ар соответствует отрицательное 
приращение объёма жидкости &И„, и наоборот. 

Уравнение (1) позволяет определить объёмную жёсткость рабочей 
жидкости, ограниченной абсолютно жёсткой оболочкой. Представим такую 
оболочку с учетом того, что из рассматриваемого объёма (рис. 1) жидкость 
может, как выходить (жидкость расширяется), так и входить (жидкость 
сжимается), то есть 


АУ = АУ,ых 75 АУ„х ' (3) 
тогда уравнение (2) примет вид 
р = О (АУ гы ах у (4) 


где аИ.„х — объём жидкости, выходящей из рассматриваемого объёма И, за 
время 4 @И», — объём жидкости, вошедшей за то же время. 


ау 
Разделив уравнение (4) на а имея в виду, что р = О - объём- 
1 


ный расход жидкости, и переставив при этом местами величины, стоящие в 
скобках, окончательно запишем: 
р г С (О 2 Ох а: ' (5) 
где ар - приращение давления жидкости в рассматриваемом объёме И за 
время аЁ О. и О,„х — суммарные расходы входящей в рассматрива- 
емый объём и выходящей из него жидкости соответственно; 
С» - объёмная жёсткость жидкости, ограниченной объёмом И, рас- 
сматриваемого участка трубы. 





Рис.1. Схема изменения объёма упругой жидкости 
при одновременном её перемещении 


Уравнение (5) позволяет определить изменение давления жидкости 
в некотором объёме, ограниченном абсолютно жёсткой оболочкой, при не- 
стационарном течении жидкости через него. Однако при использовании 
уравнения (5) следует помнить, что оно справедливо только в тех случаях, 
когда изменение давления во всех точках рассматриваемого объёма проис- 
ходит одновременно (что вполне допустимо при скоростях значительно 
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меньших, чем скорость распространения звука в среде конкретной рабочей 
жидкости). 
Объёмная жёсткость цилиндрической стенки. В случае, если дефор- 
мацией оболочки пренебречь нельзя, то изменение её объёма (например, 
объёма цилиндрической трубы), можно определить, используя следующие 
положения. 

Приращение объёма трубы при её радиальной деформации можно 
определить по формуле: 


АУ» = 2таг, 
имея при этом в виду, что 
> 
дг=—— 40, (6) 
Е 


ст 
где ги [-— соответственно радиус и длина рассматриваемого участка тру- 
бы; Е. - модуль упругости материала стенки трубы; &/-— приращение 
радиуса трубы, соответствующее приращению напряжения сб, воз- 
никающего в её стенке. 
Известно, что напряжение, возникающее в стенке трубы при воз- 
действии на неё изнутри давления р, можно рассчитать по формуле 


Г 5 
б = р, где б — толщина стенки трубы; тогда его полный дифференциал 


| 
определяем по выражению ас = 5((ар+ ра»), но так как для металли- 


ческих труб 





тар >> ра!, (7) 
то, пренебрегая вторым членом в скобке полного дифференциала напря- 
2 
т 
жения, уравнение (6) можно записать в виде 4г=— — ар. Тогда окон- 
ст 
чательно получим @=С.АУ„,, где объёмная жёсткость стенки цилин- 
дрической трубы С. определяем по формуле 
4Е 
Ся (8) 
аз! 
б 


где @- внутренний диаметр стенки трубы в ненапряжённом состоянии. 
Определяя объёмную жёсткость цилиндрической стенки по уравне- 
нию (8), следует помнить об условии (7), из которого вытекает, что урав- 
нение (8) можно использовать только для расчёта цилиндрических труб, 
изготовленных из металла, либо из другого материала с достаточно боль- 
шим модулем упругости. 
Приведённая объёмная жёсткость п гидравлических элементов, 
включённых в систему параллельно. Рассмотрим случай, когда к дан- 
ной точке гидросистемы подсоединены п различных элементов, каждый из 
которых обладает собственной объёмной жёсткостью. Тогда, допуская, что 
изменение давления в данной точке мгновенно передаётся всем точкам 
рассматриваемого объёма, можно записать: 
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а = > ау, р (9) 


где а; — изменение объёма жидкости, соответствующее деформации его 
+го элемента. 
Тогда приведённую объёмную жёсткость системы в целом опреде- 
ляем по уравнению: 








с _- а _ 1 1 
пр п п ПМ ' 
хауУ х а 
= = фр =С,; 


или после преобразований окончательно получаем: 


Пс 
СЕ, (10) 


>| Ис 


ни \ Сре 
где С„- приведённая объёмная жёсткость гидравлической системы, состо- 
ящей из п элементов, включённых в систему параллельно; 
С; — объёмная жёсткость /-го элемента системы. 
Приведённая объёмная жёсткость трубы с рабочей жидкостью. На 
практике деформация рабочей жидкости происходит одновременно с де- 
формацией стенок трубы, по которой она движется, но тогда в соответ- 
ствии с уравнением (10) можно записать 
© ИВ Оба 
о 
ж ст 
После преобразований уравнения (11) с учётом (1) и (8) получим 
выражение для определения приведённой объёмной жёсткости цилиндри- 
ческой трубы с рабочей жидкостью 


С 


(11) 


Е 

пр.тр = Е Е ь (12) 

У В 

ст 

Уравнение (12) можно также получить, подставив в (1) значение 
приведённого модуля упругости трубы [3, с.246], имея при этом в виду, что 
и = Ир . 
Приведённая объёмная жёсткость поршневой и штоковой поло- 
стей гидравлического цилиндра. Аналогично выводу уравнения (12) 
легко доказать, что приведённая объёмная жёсткость поршневой полости 
гидравлического цилиндра может быть определена по формуле: 





4Е 
Е х } (13) 


л4?х Во + 4У р 
б 


ст 





где а@- внутренний диаметр гильзы гидроцилиндра в ненапряжённом со- 
стоянии; х- ход поршня гидроцилиндра; И, „р. —- объём «вредного» 
пространства поршневой полости гидроцилиндра. 
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Аналогично приведённая объёмная жёсткость штоковой полости 
гидроцилиндра может быть определена по уравнению 


4Е (14) 
З 7 
л(1-х аа я ТА 


ст 





С при = 


где / — полный ход поршня гидроцилиндра; а„- диаметр штока гидроци- 
линдра; И, „р — объём вредного пространства штоковой полости гид- 
роцилиндра. 
Приведённая объёмная жёсткость рабочей полости плунжерного 
гидравлического цилиндра. Приведённая объёмная жёсткость рабочей 
полости плунжерного гидроцилиндра может быть определена с использо- 
ванием формул (1), (8) и (10). 
При этом имеем в виду, что для плунжерного гидроцилиндра 
Ур =У, -Ущ , (15) 


где И, и И) - соответственно внутренний объём гильзы цилиндра (с учё- 
том её вредного объёма) и объём части плунжера, находящейся в 
гильзе. 
Тогда после проведения соответствующих преобразований легко 
получить формулу для определения приведённой объёмной жёсткости ра- 
бочей полости плунжерного гидроцилиндра С;р.лл в виде 








Е 
Се = аЕ ^ ' (16) 
У — —^ — Аы(Г-х 
еее) 
ИЛИ 
4Ех 
а = Ре Е ' (17) 
п Га. (1-х) |+ 4У д 
Е, 


где (и 4, - соответственно внутренний диаметр гильзы и диаметр плун- 
жера гидроцилиндра в ненапряжённом состоянии; Л» — площадь попереч- 
ного сечения плунжера гидроцилиндра; Илл.вр - объём вредного простран- 
ства рабочей полости плунжерного гидроцилиндра. 

Заключение. Применение уравнений (1), (12)-(14), (16) и (17) совместно 
с уравнением (2) позволяет производить динамический расчёт ГМС не 
только с простыми, но и с разветвлёнными гидравлическими системами [6, 7]. 
При этом следует отметить, что объёмные жёсткости гидравлических ци- 
линдров не постоянны, так как они зависят от перемещения, скорость из- 
менения которого, в свою очередь, зависит от давления. Это делает урав- 
нения динамики ГМС нелинейными, а потому для их решения целесообраз- 
но использовать численные методы решения дифференциальных уравне- 
ний. 

На рис.2, 3 приведены графики изменения скорости поршней двух 
гидроцилиндров с диаметрами поршней 100 мм, штоков 40 мм и ходом 
поршней 0,9 м, синхронизируемых посредством дроссельного делителя пото- 
ков [7], при расходе рабочей жидкости на входе делителя 0=0,8*103 мз/с. 
Модуль упругости принимался для рабочей жидкости равным 1200 МПа, 
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для стенок трубопроводов и гидроцилиндров — 5000 МПа. Толщина стенок 
трубопроводов принималась равной 1 мм, а гидроцилиндров - 5 мм. К 
поршням гидроцилиндров приложены постоянные силы А, = АР =50 кН, 
массы перемещаемых элементов системы М =500 кг, № =50 кг. Кроме то- 
го, на поршни обоих гидроцилиндров воздействует гармоническое возбуж- 
дающее усилие А, =5 $т(8 КН, где & = 50 с"; {- время. 





Рис.2. Изменение скорости движения поршня первого гидроцилиндра 
во времени с учётом его объёмной жёсткости (1) и без учёта (2) 


Расчет производился с учётом приведённых объёмных жёсткостей 
гидроцилиндров (кривые 1) и без учёта приведённых объёмных жёсткостей 
гидроцилиндров - жёсткость принималась бесконечной (кривые 2). Объём- 
ная жёсткость прочих элементов гидросистемы (трубопроводов) учитыва- 
лась в обоих случаях. 





Рис.3. Изменение скорости движения поршня второго гидроцилиндра 
во времени с учётом его объёмной жёсткости (1) и без учёта (2) 
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Анализ изменения скорости во времени (см.рис.2,3) указывает на 
существенное отличие характера перемещения поршней гидроцилиндров 
при учёте и без учёта приведённых объёмных жёсткостей гидроцилиндров. 
Главным отличием является то, что при учёте приведённой объёмной 
жёсткости гидроцилиндров по мере выдвижения поршней амплитуда коле- 
бания их скорости увеличивается: это объясняется уменьшением приве- 
дённой объёмной жёсткости гидроцилиндров. В то же время при бесконеч- 
ной объёмной жёсткости гидроцилиндров амплитуда колебания скорости 
поршней остаётся неизменной. 

При изменении приведённой объёмной жёсткости ГМС изменяется и 
собственная частота её колебаний, что также частично объясняет рост ам- 
плитуды изменения скоростей поршней гидроцилиндров. Но это наводит на 
мысль о возможности гидромеханического резонанса, вызываемого совпа- 
дением собственноё частоты системы, при изменении приведённой объём- 
ной жёсткости её гидравлического привода, с частотой вынужденных коле- 
баний, а следовательно, заслуживает специальных исследований. 
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